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The critical behaviour of three-dimensional semi-infinite 
Ising ferromagnets at planar surfaces with (i) random surface- 
bond disorder or (ii) a terrace of monatomic height and 
macroscopic size is considered. The Griffiths-Kelly-Sherman 
correlation inequalities are shown to impose constraints on 
the order-parameter density at the surface, which yield upper 
and lower bounds for the surface critical exponent 0i. If the 
surface bonds do not exceed the threshold for supercritical en- 
hancement of the pure system, these bounds force /3i to take 
the value f3° Td of the latter system's ordinary transition. This 
explains the robustness of /3i rd to such surface imperfections 
observed in recent Monte Carlo simulations. 
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In a recent paper Pleimling and Selke (PS) Q] re- 
ported the results of a detailed Monte Carlo analysis 
of the effects of two types of surface imperfections on 
the surface critical behaviour of d = 3 dimensional semi- 
infinite Ising models with planar surfaces and ferromag- 
netic nearest-neighbour (NN) interactions: (i) random 
surface-bond disorder and (ii) a terrace of monatomic 
height and macroscopic size on the surface. For type 
(i), both the ordinary and special transitions were stud- 
ied. They found that the asymptotic temperature depen- 
dence of the disorder-averaged surface magnetization on 
approaching the bulk critical temperature T c from be- 
low could be represented by a power law ~ \t\^ 1 with 
t = (T — T c )/T c , where j3\ agreed, within the available 
numerical accuracy, with the respective values (3° rd ~ 0.8 
and /J/ 1 ~ 0.2 of the pure system's ordinary and spe- 
cial transitions. For type (ii), where the interaction con- 
stants were chosen such that only an ordinary transition 
could occur, the same value (3° rA of the perfect system 
was found for j3\. 

Their findings for the case of (i) are in conformity with 
the relevance/irrelevance criteria of Diehl and Niisser 
P,pl according to which the pure system's surface critical 
behaviour should be expected to be stable or unstable 
with respect to short-range correlated random surface- 
bond disorder depending on whether the surface specific 
heat C\\ Q of the pure system remains finite or diverges 
at the transition. It is fairly well established H|| that 
Cn approaches a finite constant at the ordinary transi- 
tion, but has a leading thermal singularity ~ | T |( d - 1 ) 1 '- 2 * 
at the special transition, where $ is the surface crossover 
exponent. In the latter case, the condition for irrelevance, 



$ < (d - 1)^/2, reduces to 

$<z/ (1) 

in d = 3 bulk dimensions. Since various Monte Carlo 
simulations |7]-[| (though not all [jl0|) and renewed field- 
theory estimates [|llj suggest a value of $ between 0.5 and 
0.6, definitely smaller than the accepted value 0.63 of v 
for d — 3, one may be quite confident that the condition 
([j]) holds. Thus short-range correlated surface-bond dis- 
order should be irrelevant in the renormalization-group 
sense at both transitions. 

Irrelevance criteria of the above Harris type (||| seem 
to work quite well in practice. Yet, from a mathemati- 
cal point of view, they are rather weak because they are 
nothing but a necessary (though not sufficient) condition 
for stability of the pure system's critical behaviour. 

In this note, I shall employ the Griffiths-Kelly-Sherman 
(GKS) inequalities Q to obtain upper and lower bounds 
on the surface magnetization densities of both types of 
imperfect systems, bounds that are given by surface mag- 
netizations of analogous systems without such imperfec- 
tions. Their known asymptotic temperature dependence 
near T c will then be exploited to obtain restrictions on 
the surface critical behaviour of the imperfect systems 
considered. For some cases of interest studied by PS [Q, 
the equality f3i = (3° rd will be rigorously established. 

Following these authors, let us consider an Ising model 
with ferromagnetic NN interactions on a simple cubic lat- 
tice of size L x x L y x L z . Periodic boundary conditions 
will be chosen along two principal axes (the x and y di- 
rections), and free boundary conditions along the third 
one (the z direction), so that the surface consists of the 
top layer at z = 1 and the bottom layer at z = L z . As- 
sociated with each pair of spins on NN sites i and j is 
an interaction constant J(i,j) > 0, which we assume to 
have the same value J whenever i or j (or both) belong 
to layers with 1 < z < L z . 

In the case of surface-bond disorder, which we consider 
first, the J(i,j) = J (s) (i,j) of all NN pairs of surface 
sites are independent random variables. The probability 
density P(Ji) of any one of these will be assumed to have 
support only in the interval [J± , J{~] (with J{" > J< > 
0). This is in conformity with, but less restrictive than, 
PS's assumption that J\ takes just two values and 
J^, either one with probability 1/2. We will also assume 
that all (bulk and surface) spins are exposed to the same 



f 



magnetic field H > 0, whose limit H — > + will be taken 
after the thermodynamic limit has been performed. 

Let K = J/k B T and h = H/J. Define to be 
the set of all dimensionless surface coupling constants 
J (h i) I J ■ Let m(i; K, , h) = be the thermal av- 
erage of a spin at site i for a given disorder configuration 
r( s \ and denote the corresponding quantity of the per- 
fect system with uniform NN surface coupling J± = rJ as 
m(i; K, r, h). Since all interactions are ferromagnetic, the 
GKS inequalities |l2| are valid. Averages of products of 
spin variables are monotone non-decreasing functions of 
all variables J(i,j) and H. Hence, for finite L x , L y , and 
L z , m(i; K, r^ s \ h) is bounded by m(i; K, r < , h) from be- 
low and by m(i; K, r > , h) from above. We choose i = i s 
to be a surface site, take the thermodynamic limit (first) 
and then let H — > + . The bounds converge towards the 
respective values of mi(K, r, + ), the spontaneous mag- 



netization of the surface layers per site, for r 
r > . Thus we obtain 



and 



TOi(if,r < ,0+) < m(i s ;K,r is \0 + ) < m 1 (K,r > ,0 + ) . 

(2) 

The following limiting forms of mi are well established 



dlrf^ll + o(t)} asr 
'"j = C^T^fl + o(r)j as r 
mic -I- 0{t) as r 



0~ at fixed r < r c 
0~ at fixed r = r c 
0^ at fixed r > r r 



(3) 



where r c ~ 1.50 is the critical value associated with 
the special transition. The quantities m\ c > 0, C%, and 
C{ are nonuniversal, whence the first two depend on r. 

Consider first the case r> < r c . Let C < and C > be 
the values of C\ for r — r < and r = r > , respectively. 
(These satisfy < C < < C> < oo provided < J < oo 
and < < J i < oo.) It follows that there exists a 
number e > independent of the disorder configuration 
r( s ) such that 



C> < m[i s ;K(T),r( s \0 + ] It^^ < C ; 



(4) 



whenever — e < r < 0. We denote the average of a quan- 
tity Q over all choices of the random variables as Q. 
Upon averaging m(i s ; .) to obtain the disorder-averaged 
surface magnetization TnT, we see that the inequality (^) 
holds for rrTf |r| - ' 31 as well. An elementary consequence 
is: If mT has a well-defined critical exponent f3f ls in the 
sense that Ml 



dis 



in 



exists, then we have 



lim 



lnmi(r) 
lnr 



ord 



(5) 



(6) 



Two further implications of (Q) are worth mentioning. 
First, if a surface critical exponent f3f ls can be defined via 
the analog of (||) for the most probable value of m(i s ; .) 
fuH , then it must have the same value /3° Td . Second, the 
inequality (Q) also rules out a limiting t dependence of 
the form ~ \t\P~ 1 \ In |t| \ v (standard logarithmic correc- 
tions) for TfTi and the most probable value of m(i s ; .). 

Consider next the case r > — r c . Let us again make 
the assumption that the limit (||) or the analogous one 
defining [3f ls exist. Then the inequalities 



PT < Pi s < PT A 



(7) 



and their analogs for j3f la can be deduced from (^). (Cf. 
Lemma 3 of fl5j.) 

The same reasoning applied in the case r > > r c shows 
that Pf ls or Pf ls must obey the relations 



< /3f s < 



(8) 



whenever the limits (^|) through which we defined them 
exist. 

Likewise in the case r < — r c , the possible values of (3f ls 
or ftf ls are restricted by 



< pf is < pl p 



(9) 



at transitions at which TrTi or the most probable value 
of m{i s ; .) Jl6[ approach zero, respectively. On the other 
hand, it should be recalled that the surface critical ex- 
ponent /?J X of the pure system's extraordinary transition 
requires a definition other than (||) : One must subtract a 
regular background contribution m^ cs from mi and define 
/3® x through the limiting behaviour mi — m^ cg ~ {t^ 1 . 
For transitions of the impure systems at which mi ap- 
proaches a constant ^ 0, it would also not make much 
sense to define /3f m via (|3|). Of course, for surface criti- 
cal exponents /3f ls not given by (||) , the above bounds do 
not apply. This means that they cannot be utilized to 
draw conclusions about the surface critical exponent /3^ x 
of the impure system's extraordinary transition. How- 
ever, for a special transition of the impure system with 
mT(r = 0) = 0, the inequalities (Js[) hold. 

The inequality (0) rules out that the impure system has 
an ordered surface phase for T > T c whenever r > < r c . 
In order that the impure system can have an extraordi- 
nary or special transition, the distribution P{J\) of the 
surface couplings typically will have to extend beyond the 
critical-enhancement threshold r c J of the pure system. 
But even if r > > r c , an ordinary transition may still oc- 
cur if the surface bonds 'on average' are not sufficiently 
enhanced (cf. jjj). However, if P{J\) extends beyond 
r c J, then disorder configurations for which macroscop- 
ically large surface regions have the same supercritical 
value (> r c J) of J\ occur with finite probability. This 
happens even if the impure system (for a typical realiza- 
tion of disorder) undergoes an ordinary transition, albeit 



2 



with exponentially small probability. By analogy with 
the bulk case fL7| , I expect surface quantities like W\ 
and the disorder-averaged surface free energy to be non- 
analytic functions of the surface magnetic field Hi at 
Hi = for temperatures between the bulk critical tem- 
perature T c and the temperature T s (' ,> ) > T c at which 
the semi-infinite pure system with homogeneous surface 



coupling Ji 



'J undergoes a transition to a surface- 



ordered, bulk-disordered phase. That is, they should dis- 
play Griffiths singularities fll7| , a problem on which we 
will not embark further here. 

Turning now to the case of surfaces with a terrace, 
we start from a pure Ising model of the sort considered 
above. Just as PS, we assume that all NN couplings 
J(i,j) (including those between surface sites) have the 
same value J. Let us denote thermal averages pertain- 
ing to this system by a superscript [I], writing, e.g., 
m^(i;K) — (si)^. We consider another system, [II], 
which differs from [7] through the addition of a zeroth 
layer at z = whose spins are assumed to interact among 
themselves and with the spins in the z = 1 layer via NN 
interaction constants Ji and J, respectively. To obtain 
a system with a terrace, [T], we choose a subregion of 
the zeroth layer (the terrace) and remove all those NN 
bonds J and Ji that are connected to lattice sites of this 
layer outside the terrace region. PS considered a strip- 
like terrace of size (L x /2) x L y , and assumed that Ji = J. 
For our considerations, the precise form and size of the 
terrace region will not be important. (One could even as- 
sume that an arbitrary subset of the spins in the zeroth 
layer are decoupled from the rest of the system.) 

Let %i be an arbitrary lattice site in the z = 1 layer. 
Since the systems [I] , [T] , and [II] differ by the addition 
of ferromagnetic interactions, we have from the GKS in- 
equalities, 



m 



where, as before, 
netic field and r = 



h = H/J > is a uniform mag- 
= Ji/J. In the thermodynamic limit 



L x , I y , L z — > oo, the lower and upper bounds converge 
towards mi(K, h), the magnetization per site of the top- 
most layer, and to mi(K, r, h), the magnetization per site 
of the layer underneath the topmost layer, respectively. 
If we assume that r < r c (subcritical surface enhance- 
ment) and take the limit h — > + , then the limiting form 
shown in the first line of (^) applies to both mi and 7712 
(with different values of Ci). As a straightforward con- 
sequence we find that the surface critical exponent fii 
of mP^(i\\ K, r, + ) (for an arbitrary site i\ with z = 0) 
strictly satisfies f3\ = /3° rd . 

It evident that the same reasoning can be applied 
to the analogous two-dimensional model with a ter- 
race to conclude that (3i takes the exactly known value 
/3° ld — 1/2. Likewise, the inequality (|J) and the re- 
sult carry over to the two-dimensional case, giving 



f3f ls = 1/2 for all values of r < 00, since r c = 00 for 
d = 2. Note also that the inequality (Q) excludes the 
possibility of an asymptotic temperature dependence of 
the form rh~i » const | -7- 1 1 / 2 1 In |r|| p (i.e., of logarithmic 
correction factors). This is because it is known for the 
pure case that no such logarithmic corrections appear in 
the limiting form of mi. 

Results of Monte Carlo simulations on the surface criti- 
cal behaviour of two-dimensional Ising models with bond 
disorder have been reported in two recent papers [ fLs] . 
However, in this work random bond disorder was as- 
sumed to be present both in the bulk and at the surface, 
a case not captured by our reasoning. Nevertheless, mi 
was found to behave as i") 1 / 2 , apparently without loga- 
rithmic corrections, even though the presence of such a 
correction could be detected in the limiting form of the 
disorder-averaged bulk order parameter. 

I am indebted to W. Selke for informing me about 
the work jl| prior to publication, and to him, Joachim 
Krug, and Kay Wiese for a critical reading of the 
manuscript. This work has been supported by the 
Deutsche Forschungsgemeinschaft through the Leibniz 
program. 
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